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We study the properties of the Google matrix of an Ulam network generated by intermittency 
maps. This network is created by the Ulam method which gives a matrix approximant for the Perron- 
Frobenius operator of dynamical map. The spectral properties of eigenvalues and eigenvectors of 
this matrix are analyzed. We show that the PageRank of the system is characterized by a power 
law decay with the exponent j3 dependent on map parameters and the Google damping factor a. 
Under certain conditions the PageRank is completely delocalized so that the Google search in such 
a situation becomes inefficient. 

PACS numbers: 05.45.-a, 89.20.Hh, 05.45.Ac 



I INTRODUCTION 

In 60s Ulam proposed a method to construct a matrix 
approximant for a Perron-Frobenius operator of dynami- 
cal systems which is now known as the Ulam method [1] . 
The Ulam conjecture was that, in the limit of small cell 
discretization of the phase space, this method converges 
and gives the correct description of the Perron-Frobenius 
operator of a system with continuous phase space. This 
conjecture was shown to be true for hyperbolic maps of 
the interval [2] . Various types of more generic maps of an 
interval were studied in [3, 4, 5]. Further mathematical 
results have been obtained in [6, 7, 8, 9] with extensions 
and prove of convergence for hyperbolic maps in higher 
dimensions. The mathematical analysis of non- uniformly 
expanding maps is now in progress [11] . At the same time 
it is known that the Ulam method applied to Hamilto- 
nian systems with integrable islands of motion destroys 
the invariant curves thus producing a strong modification 
of properties of the Perron-Frobenius operator of the sys- 
tem with continuous phase space (see e.g. [12]). 

Recently it was shown that the Ulam method natu- 
rally generates a class of directed networks, named Ulam 
networks, which properties have certain similarities with 
the World Wide Web (WWW) networks [12]. Thus the 
Google matrix constructed for the Ulam networks built 
for the Chirikov typical map has a number of interesting 
properties showing a power law decay of the PageRank 
vector. 

The classification of network nodes by the PageRank 
Algorithm (PRA) was proposed by Brin and Page in 1998 
[13] and became the core of the Google search engine used 
everyday by majority of internet users. The PRA is based 
on the construction of the Google matrix which can be 
written as (see e.g. [14] for details): 

G = aS + (1 - a)E/JV . (1) 

Here the matrix S is constructed from the adjacency ma- 
trix A of directed network links between N nodes so that 
Sij = Aij I ^2 k Akj and the elements of columns with only 



zero elements are replaced by 1/N. The second term in 
r.h.s. of (1) describes a finite probability 1 — a for WWW 
surfer to jump at random to any node so that the matrix 
elements E^j — 1. This term stabilizes the convergence 
of PRA introducing a gap between the maximal eigen- 
value A = 1 and other eigenvalues A^ . Usually the Google 
search uses the value a — 0.85 [14]. The factor a is also 
called the Google damping factor. By the construction 
Si @ij — 1 so that the asymmetric matrix G belongs to 
the class of Perron-Frobenius operators. Such operators 
naturally appear in the ergodic theory [15] and dynam- 
ical systems with Hamiltonian or dissipative dynamics 
[16, 17]. 

The right eigenvector at A = 1 is the PageRank vector 
with positive elements pj and J* 1 , pj = 1, the components 
Pj of this vector are used for ordering and classification 
of nodes. The PageRank can be efficiently obtained by a 
multiplication of a random vector by G which is of low 
cost since in average there are only about ten nonzero 
elements in a typical line of G of WWW. This procedure 
converges rapidly to the PageRank. All WWW nodes 
can be ordered by decreasing pj > pj+i) so that 
the PageRank plays a significant role in the ordering of 
websites and information retrieval. The classification of 
nodes in the decreasing order of pj values is used to clas- 
sify importance of network nodes as it is described in 
more detail in [14]. 

Due to a spectacular success of the Google search the 
studies of PageRank properties became very active re- 
search filed in the computer science community. A num- 
ber of interesting results in this field can be find in 
[18, 19, 20, 21]. An overview of the field is available 
in [22] . It is established that for large WWW subsets pj 
is satisfactory described by a scale-free algebraic decay 
with pj ~ 1 / j@ where j is the PageRank ordering index 
and (3 « 0.9 [14, 23]. 

In this work we analyze the properties of Google ma- 
trix constructed from Ulam networks generated by one- 
dimensional (ID) intermittency maps. Such maps were 
introduced in [24] and studied in dynamical systems with 
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FIG. 1: (Color online) Two types of intermittency map of the 
interval given by the map functions j\(x) and ji{x), the func- 
tions are identical for < x < 1/2 but have different branches 
at 1/2 < x < 1 with fi(x) (red/gray) and f2(x) (blue/black); 
map functions are shown at z\ = 2 and at different values of 
parameters 22 and a; the straight line shows f(x) — x. 



intermittency properties (see e.g. [25, 26, 27, 28]). A 
number of mathematical results on the measure distribu- 
tion and slow mixing in such maps can be find in [29, 30] 
(see also Refs. therein). The mathematical properties of 
convergence of the Ulam method in such intermittency 
maps are discussed in a recent work [11]. The analysis 
of such ID maps is simpler compared to the 2D map 
considered in [12]: for example the PageRank at a = 1 
is described by the invariant measure of the map which 
can be find analytically as a function of map parameters. 
Following the approach discussed in [12, 31] we study 
not only the PageRank but also the spectrum and the 
eigenstates of the Google matrix generated by the inter- 
mittency maps. Indeed, the right eigenvectors tpi and 
eigenvalues Xi of the Google matrix (Gtpi — A^i) are 
generally complex and their properties should be studied 
in detail to understand the behavior of the PageRank. 
We show that under certain conditions the properties of 
the PageRank can be drastically changed by parameter 
variation. 

The results are presented in a following way: in Sec- 
tion II we describe the class of intermittency maps and 
the distribution of links in the corresponding Ulam net- 
work, the spectral properties of the Google matrix and 
PageRank are considered in Sections III and IV, the dis- 
cussion of the results is presented in Section IV. 



FIG. 2: (Color online) Google matrix at a = 1 generated by 
the intermittency map fi(x) at 21 = 2, 22 = 0.2, N = 50, 
N c — 10 6 (amplitude of matrix elements is changing from 
zero (black/blue) to 1 (red/gray)). 



INTERMITTENCY MAPS 

The intermittency maps of the interval considered in 
this paper are described by the two map functions de- 
pending on parameters and defined for the first model 
as: 



x + (2x) Zl /2, for < x < 1/2 
h{x) = { (2x-l-(l- x) Z2 + 1/2 Z2 )/(1 + 1/2 22 ) 
for 1/2<x<1 



and for the second model as 



/a(z) 



x+ (2x) Zl /2, for Q<x<l/2 
asin[7r(x-l/2)], for 1/2 < x < 1 



(2) 



(3) 



The parameters Zi,Z2,a are positive numbers. The dy- 
namics is given by the map x — fi{x) and x = /^(x)- 
The map functions fi t i(x) are shown in Fig. 1. 

According to the usual theory of intermittency maps 
and ergodicity theory [24, 25, 26, 29, 30, 32, 33] in the 
case of chaotic dynamics the steady state invariant dis- 
tribution g(x) of the map is proportional to a time t(x) 
spent by a trajectory at point x which is proportional to 
t ~ l/a; 1_Zl so that one has a power law distribution at 
small values of x: 



g(x) oc l/x zi ~ 



(4) 



F° r /i-map the dynamics is fully chaotic while for fa- 
map a fixed point attractor appears for a > 0.945 when 
f 2 (x) = x. 
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FIG. 3: (Color online) Distribution of links between nodes 
of the Ulam network of N = 10 6 size for the first model 
with z\ = 2, Z2 — 0.2 (top panel) and the second model with 
z\ = 2, a = 0.9 (bottom panel). Here Nl(k) gives the num- 
ber of nodes which have k outgoing (black points) or ingoing 
(red/gray squares) links respectively. Insets show data for 
small k values in linear scale. The straight line shows the the- 
oretical slope for outgoing links (Nl(k) oc 1/k 9//4 first model, 
top panel) and for ingoing links (Nl(k) oc 1/k 3 second model, 
bottom panel). 



The Ulam networks generated by the intcrmittency 
maps (2), (3) are constructed in a way similar to one 
described in [1, 12]: the whole interval < x < 1 is 
divided on N equal cells and N c trajectories (randomly 
distributed inside cell) are iterated on one map iteration 
from cell j to obtain matrix elements for transitions to 
cell i: Sij = Ni(j)/N c where Ni(j) is a number of tra- 
jectories arrived from cell j to cell i. The image of the 
density of Google matrix elements is shown in Fig. 2 for 
the first model. The structure of the matrix repeats the 
form of the map function fi(x). We used from 10 4 to 10 6 
cell trajectories N c , the obtained results are not sensitive 
to N c variation in this interval. 

The differential distribution of number of nodes Nl(k) 
with ingoing or outgoing links n is shown in Fig. 3. The 
first model shows a sharp drop of ingoing links and a 
power law decay of outgoing links. For the second model 
the situation is inverted. These properties can be under- 
stood from the following arguments. For the first model, 
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FIG. 4: (Color online) Distribution of eigenvalues A in the 
complex plain for the Google matrix at a = 1 for the first 
(zi — 2, Z2 = 0.2, left panel) and second (zi = 2, a = 0.9, 
right panel) models at N = 12000. Color of small squares is 
determined by the value of PAR £ associated with the corre- 
sponding eigenvector ipi as show in the palette (the values of 
£ are averaged over the states inside of the square size). 



the number of outgoing links is k = 
the derivative is diverging near x 



dx/dx = dfi{x)/dx, 
— 1 where we have 



x 



is N n ~(l- x) 
bution of nodes 



(i z 2 )^ rp ne num ber of nodes with k links 
I/k 1 ^ 1-22 ) and the differential distri- 



N' 



L ~ dN n /dx ~ l//s", n = (2 - z 2 )/(l - z 2 ) 



(5) 



For the data of Fig. 3 (top panel) at z 2 = 0.2 this estimate 
gives /x = 9/4 in a good agreement with the numerical 
data. For the second model df2(x)/dx is always finite 
and we have a sharp drop for outgoing links distribution. 
The number of ingoing links is n = dx/dx ~ l/x 1 ^ 1 ^ 2l/ 
since we have x ~ (1 — x) 2v near x = 1 (in our case v = 1 
but we consider here a general case). Hence, the number 
of nodes with k links is N n ~ x ~ 1/ ' n 2u /^ 2 "^ 1 ) and 



dN n /dK~l/K», n= (4i/-1)/(2i/-1) . (6) 



For our case with v = 1 we have fi = 3. This value is in 
a good agreement with the data of Fig. 3. For the first 
model dx/dx is always finite and we have a sharp drop 
of ingoing links distribution. 

This analysis allows to understand the origin of power 
law distributions of links in the Ulam networks generated 
by Id maps. 



SPECTRAL PROPERTIES OF THE GOOGLE 
MATRIX 

The distribution of the eigenvalues of the Google ma- 
trix at a = 1 constructed from the Ulam network de- 
scribed above is shown in Fig. 4 for two models (2) 
and (3). As in [12, 31] we characterize an eigenstate 
ipi by a PArticipation Ratio (PAR) defined as £j = 
(Ej Mi)\ 2 ) 2 /Ej In fact PAR gives an effec- 

tive number of nodes populated by a given eigenstate, it 
is broadly used in systems with disorder and Anderson lo- 
calization. The states are normalized by the condi- 
tion J^j l^»(i)| 2 = !• F° r the PageRank pj proportional 
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FIG. 5: (Color online) Dependence of density of states W(-y) 
on 7 shown for different values of N for the first (at z\ — 
2, Z2 = 0.2, top panel) and second (at z% = 2, a = 0.9, bottom 
panel) models; G matrix is taken at a = 1. 



to -01 (j), ordered in the decreasing order of probability, 
we use also probability normalization ^2jPj = 1- 

There are few main features of the spectrum of A in 
Fig. 4 visible for two models: there are states with |A| 
close to 1 which have relatively small values of £; there 
is a circle like structure of eigenvalues and the maximum 
PAR are in the middle ring around the center. The large 
circle is present for both maps fi{x) and f2(x). This 
means that it appears due to the left branch of the map 
corresponding to intermittent motion near x = 0. The 
density distributions W(j) = dN^/dj in the decay rate 
denned as 7 = — 2 In | A| are shown in Fig. 5 (here dN 1 is 
a number of states in the interval d'y). It is clear that in 
the limit of large matrix size N we have a convergence 
to a limiting distribution which has a characteristic peak 
at 7 w 2. 

Examples of few eigenstates ip(i) with values of 7 m = 
—2 In |A m | equal and close to zero are shown in Fig. 6 (the 
index 1 < i < N gives the cell position Xi = (i — 1)/N, 
index m orders j m from zero to maximum 7). The first 
state with Ai = 1 is the steady state distribution 

generated by the map /1 (x) (the states for the map ji (x) 
have similar structure and we do not show them here). 
We have oc 1/i 13 with [3 = 1 for Z\ = 2 is agreement 
with the theoretical expression (4) (the numerical fit gives 
/3 = 0.97). The state ipi(i) is monotonic in i so that it 
coincides with the PageRank pj up to a constant factor. 
Eigenstates with next values of 7 are characterized by the 
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FIG. 6: (Color online) Absolute value of three eigenstates ip(i) 
for the first model /i(ar) with zi = 2, z 2 = 0.2 and N = 12000. 
The eigenstates correspond to the eigenvalues Ai = 1.0 (black 
circles), A2 ~ 0.9998 (red squares) and A4 ~ 0.9983 (green 
diamonds) (see Fig. 4, left panel). The corresponding PAR 
values are £1 ~ 2.54, £2 ~ 1-21 and £4 « 9.00 respectively. 
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FIG. 7: Same as in Fig. 6 for the eigenstates with A51 ~ 0.78+ 
i0.42 (a state in a large circle, top panel), A6i « 0.56 + i0.68 (a 
state in a large circle, middle panel) and A1010 ~ — 0.36 + i0.35 
(a state in the dense part of the spectrum, bottom panel); the 
corresponding PARs are £ ~ 1231, 1482, 4367 respectively. 



same decay at large i with additional minima at certain 
values of i similar to few nodes of eigenstates in quantum 
mechanics. 

The structure of eigenstates is changed when the value 
of 7 is increased. Typical states are shown in Fig. 7. 
The states on the first circle of |A| have peaked structure 
at certain i with a plateau at large i. For 7 values at 
the maximum of W(j) (see Fig. 5) the eigenstates are 
delocalized over the whole interval of 1 < i < N. 

An effective number of sites contributing to an eigen- 
state can be characterized by the PAR £. For the PageR- 
ank the value of £ is independent of the matrix size N 
as it is clearly shown in Fig. 8. This is due to the power 
law decay of the PageRank pj ~ 1/j which corresponds 
to an algebraic localization. The dependence of £ on 7 is 
shown in Fig. 9. For small 7 it can be fitted by a power 
law growth £ ~ 7 12 . The origin of the exponent of this 
growth requires further analysis. 

Finally we note that we also determined the depen- 
dence of number of states iV 7 with values of 7 > 5 on the 
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FIG. 8: Dependence of PAR £ of the PageRank on matrix 
size N: the circles (black) are for the first model fi(x) with 
Z\ = 2, Z2 = 0.2 and the squares (red/grey) are for the second 
model f2(x) with z\ = 2, a = 0.9. 
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FIG. 10: Dependence of gap A 12 = 1 — | Aa | between the first 
eigenstate with Ai = 1 and next one with maximum | A3 1 on 
N for the first fi(x) (21 = 2, 22 = 0.2) and second f2(x) 
(zi — 2, a = 0.9) models. The straight dashed line shows the 
dependence A12 oc 1/N. 




FIG. 9: Dependence of PAR £ on 7 for the first model with 



Zl — 2, Z2 



0.2 and N = 10 . Inset show data in log-log 



scale with growth of £ ~ 7 ' at small values of 7. 




matrix size N. Our data (not shown) are well described 
by the dependence N~ ( ~ N so that in contract to the 
results presented in [12] there are no singes of the fractal 
Weyl law. We attribute this to the fact that in contrast 
to the dissipative map with a global contraction stud- 
ied in [12] in the intermittency maps all dynamics takes 
place on the whole one-dimensional interval with inhomo- 
geneous distribution of measure but without fractality. 

PROPERTIES OF THE PAGERANK 

The spectral gap between Ai — 1 equilibrium state 
and the next state with maximum | A2 1 has very small 
gap A 12 = 1 — I Aa] which goes to zero with the increase 
of N like A 12 w 3/N (see Fig. 10). This happens due to 
the dynamical properties of the maps (2), (3) where the 
time spent at small x ~ 1/N is of the order t x ~ l/a: 21-1 
(see e.g. [24, 32, 33]), so that the corresponding A 12 ~ 
l/t x ~ l/N 21 ^ 1 that gives the exponent 1 for z\ = 2. 

Due to such decrease of A12 with TV the PRA has 
bad convergence at a = 1 for large values of AT. Up 
to N ~ 14000 we use the direct diagonalization of G 
matrix which gives an algebraic decay Pj ~ 1/j with 
[3=1 (see Fig. 6). For larger value of N we used the 
continuous map obtaining pj from an equilibrium distri- 



FIG. 11: (Color online) Dependence of the PageRank pj on 
j for the first model at different values of a (21 = 2, 22 =0.2, 
top panel) and different values of zi (22 = 0.2, a = 1, bottom 
panel). The data are obtained from the continuous map (see 
text) for a — 1 and the PageRank algorithm at a < 1, the 
number of nodes is N — 10° . 

bution over the cells of size 1 /N after a larger number of 
map iterations U « 10 9 and large number of trajectories 
Ntr ~ 10. This distribution converges to a limiting one 
at large values of U (see Fig. 11). Both methods give the 
same result for N < 2 ■ 10 4 . The numerical data for the 
exponent (3 are in good agreement with the theoretical 
dependence (4) (3 = Z\ — 1 as it is shown in Fig. 12 (we 
attribute small deviations from the theoretical values to 
finite size effects of N). 

For a < 1 the PRA, described in the Introduction, is 
stable and converges rapidly to the PageRank. It gives 
the same results as the exact diagonalization for N < 
2 • 10 4 . The dependence of PageRank on a is shown 
in Fig. 11 (top panel). A small decrease down to a = 
0.999 modifies pj at j < 100 making pj very flat in this 
region. For a = 0.875 the PageRank becomes completely 
delocalized over the whole system size N. 

For the second model the PageRank depends strongly 
on the value of a. For a < 0.945 when the dynamics 
is chaotic and the steady-state distribution is given by 
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FIG. 12: Dependence of the PageRank exponent /3 (pj ~ 
l/f) on zi for the first model at z 2 = 0.2 and a = 1, N = 10 5 . 
The straight dotted line shows the fit (3 = 1.042zi - 1.107. 
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FIG. 13: (Color online) Dependence of the PageRank pj on j 
for the second model at different values of a {z\ = 2, a — 0.96, 
top panel) and different values of a (zi = 2, a = 1, bottom 
panel). The data are obtained from the continuous map (see 
text) for a = 1 and the PageRank algorithm at a < 1, the 
number of nodes is N — 10 . 



Eq. (4) the properties of the PageRank are similar to those 
of the first model described above, e.g. we have j3 = 1 
being independent of a for a = l,z% = 2 (see Fig. 13, 
bottom panel). However, for a > 0.945 the map has a 
fixed point attractor and the Page Rank becomes local- 
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FIG. 14: (Color online) Dependence of PAR £ of the PageR- 
ank (£ values are shown by color) on parameters a and a for 
the second model at z\ = 2; N = 10 s ; arrow marks the region 
with higher resolution shown on the right panel. 



ized practically on one site at a = 1. In this regime with 
fixed point attractor the PageRank is very sensitive to a 
variation: at a < 1 we have pj ~ l/j' 3 with the fit values 
(3 w 0.79 at a = 0.98, (3 « 0.60 at a = 0.9, (3 « 0.42 at 
a = 0.8 and /3 « 0.32 at a = 0.7. 

The derealization of the PageRank from the fixed 
point attractor state is also clearly seen in the variation 
of PAR £(a, a) shown in Fig. 14. This shows that even 
if at a = 1 the PageRank is dominated only by one node 
a decrease of a allows to obtain weighted contribution of 
other nodes. 

We also note that in the phase of fixed point attractor 
the spectrum of eigenvalues A has globally a structure 
rather similar to one at a — 0.9 < 0.945 (see Fig. 4, 
right panel). However, the PAR values of all eigenstates 
at a > 0.945 become rather close to unity showing that 
almost all eigenstates are strongly localized in this phase. 
For example, for a — 0.96, we have almost all & in the 
range from 1 to 4 for N = 10 4 , it is interesting that about 
53% of the states have |A| < e~ 10 (for a — 0.9 this circle 
in A contains 23% of states, see Fig. 4 right panel). 



DISCUSSION 

The present studies allowed to establish a number of 
interesting properties of the Google matrix constructed 
for the Ulam network generated by intermittency maps. 
A general property of such networks is the existence of 
states with eigenvalues |A| being very close to unity. The 
PageRank of such networks at a — 1 is characterized by a 
power law decay with an exponent determined by the pa- 
rameters of the map. It is interesting to note that usually 
for WWW it is observed that the decay of the PageR- 
ank follows the decay law of ingoing links distribution 
N™(k) (see e.g. [21]). In our case the decay of PageR- 
ank is independent of N™(k) decay as it is clearly shown 
by Eqs. (5), (6) and the data of Figs. 3,11,13. In fact 
a map with singularities of both maps fi(x) and f2{x) 
(e.g. /s (x) which behaves like x + x Zl at small x, like 
(1/2 - x) zi at x < 1/2 close to 1 /2 and like (1 - x) v near 
x = 1) will have the asymptotic decay of links distribu- 
tion given by Eqs. (5), (6) but the decay of the PageRank 
will be given by (3 = z\ — 1, hence, being independent of 
the decay of links distribution. 

Our results also show that while at a close to unity the 
decay of the PageRank has the exponent (3 w 1 but at 
smaller values a rj 0.9 the PageRank becomes completely 
delocalized (see Fig. 11). In this delocalized phase the 
PAR £ grows with the system size approximately as £ oc 
N. The derealization of the PageRank can also take 
place at a = 1 due to variation of the parameters of 
the map (e.g. for z% — > 1). It is rather clear that the 
derealization of the PageRank makes the Google search 
inefficient. 

We hope that the properties of Ulam networks gen- 
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erated by simple maps will be useful for future studies 
of real directed networks including WWW. Indeed, the 
whole world will go blind if one day the Google search will 
become inefficient. The investigations of the Ulam net- 
works can help to understand the properties of directed 
networks in a better way that can help to prevent such a 
dangerous situation. 
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